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Abstract. A fc-noncrossing RNA pseudoknot structure is a graph over {1, . . . ,n} without 1- 
arcs, i.e. arcs of the form (i, i + 1) and in which there exists no fc-set of mutually intersecting 
arcs. In particular, RNA secondary structures are 2-noncrossing RNA structures. In this paper 
we prove a central and a local limit theorem for the distribution of the numbers of 3-noncrossing 
RNA structures over n nucleotides with exactly h bonds. We will build on the results of I10| 
and [1 1) . where the generating function of fc-noncrossing RNA pseudoknot structures and the 
asymptotics for its coefficients have been derived. The results of this paper explain the findings 
on the numbers of arcs of RNA secondary structures obtained by molecular folding algorithms 
and predict the distributions for fc-noncrossing RNA folding algorithms which are currently being 
developed. 



1. Introduction 



An RNA molecule consists of the primary sequence of the four nucleotides A, G, U and C together 
with the Watson-Crick (A-U, G-C) and (U-G) base pairing rules. The latter specify the pairs of 
nucleotides that can potentially form bonds. Single stranded RNA molecules form helical structures 
whose bonds satisfy the above base pairing rules and which, in many cases, determine their function. 
For instance RNA ribosomes are capable of catalytic activity, cleaving other RNA molecules. Not 
all possible bonds are realized, though. Due to bio-physical constraints and the chemistry of 
Watson-Crick base pairs there exist rather severe constraints on the bonds of an RNA molecule. 
In light of this three decades ago Waterman et.al. pioneered the concept of RNA secondary 
structures |21| I17j. being subject to the most strict combinatorial constraints. Any structure can 
be represented by drawing the primary sequence horizontally, ignoring all chemical bonds of its 
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backbone, see Fig. [TJ Then one draws all bonds, satisfying the Watson-Crick base pairing rules 
as arcs in the upper half-plane, effectively identifying structure with the set of all arcs. In this 
representation, RNA secondary structures have no 1-arcs, i.e. arcs of the form (i, i + 1) and no two 
arcs («2, J2), where i\ < j\ and 12 < 31 with the property i\ < 12 < ji < j%. In other words 

there exist no two arcs that cross in the diagram representation of the structure. It is well-known 




A G G C A A U £ U A CAGCGU 

A 




FIGURE 1. RNA secondary structures. Diagram representation (top): the primary 
sequence, AGGCAAUCUACAGCGU, is drawn horizontally and its backbone bonds 
are ignored. All bonds are drawn in the upper half-plane. Secondary structures have the 
property that no two arcs intersect and all arcs have minimum length 2. Outer planar 
graph representation (bottom). 

that there exist additional types of nucleotide interactions [T] . These bonds are called pscudoknots 
[2"5] and occur in functional RNA (RNAseP P3]), ribosomal RNA [T2] and are conserved in the 
catalytic core of group I introns. Pscudoknots appear in plant viral RNAs pseudo-knots and in in 
vitro RNA evolution [20] experiments have produced families of RNA structures with pscudoknot 
motifs, when binding HIV-1 reverse transcriptase. Important mechanisms like ribosomal frame 
shifting [3] also involve pseudoknot interactions, fc-noncrossing RNA structures introduced in [10] 



CENTRAL AND LOCAL LIMIT THEOREMS FOR RNA STRUCTURES 



3 



capture these pseudoknot bonds and generalize the concept of the RNA secondary structures in 
a natural way. In the diagram representation fc-noncrossing RNA structure has no 1-arcs and 
contains at most k — 1 mutually crossing arcs. 
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FIGURE 2 . fc-noncrossing RNA structures, (a) secondary structure (with isolated labels 
3,7,8,10), (b) planar 3-noncrossing RNA structure, 2,9 being isolated (c) the smallest 
non-planar 3-noncrossing structure 



The starting point of this paper was the experimental finding that 3-noncrossing RNA structures for 
random sequences of length 100 over the nucleotides A, G, U and C exhibited sharply concentrated 
numbers of arcs (centered at 39). It was furthermore intriguing that the numbers of arcs were 
significantly higher than those in RNA secondary structures. While it is evident that 3-noncrossing 
RNA structures have more arcs than secondary structures, the jump from 27 to 39 (for n = 100 
) with a maximum number of 50 arcs was not anticipated. Since all these quantities were via 
the generating functions for fc-noncrossing RNA structures in |10j explicitly known we could easily 
confirm that the numbers of 3-noncrossing RNA structures with exactly h arcs, S' 3 (n,h) satisfy 
indeed almost "perfectly" a Gaussian distribution with a mean of 39, see Fig. [3] We also found 
that a central limit theorem holds for RNA secondary structures with h arcs, see Figure HJ These 
observation motivated us to understand how and why these limit distributions arise, which is what 
the present paper is about. Our main results can be summarized as follows: 
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FIGURE 3. Central limit theorem and local limit theorem for 3-noncrossing RNA struc- 
tures of length n = 100 with exactly h arcs: we display the central limit theorem (left) 
for S^lOO, h), h = 1, 2, ■ ■ ■ 50 (labeled by red dots) with mean 0.39089 • 100 = 39.089 and 
variance 0.041565 • 100 = 4.1565, and for the local limit theorem (right), we display the 
difference V4.1565 P ( ~ y "~ 39 ' 089 = x ] }=e~^r which is maximal close to the peak of 

\ V4.1565 / V2vr v 

the distribution. 



Theorem. Let S' 3 (n,h) denote the number of 3-noncrossing RNA structures with exactly h arcs. 
Then the random variable X n having distribution P(X n = h) = S' 3 (n,h)/S3(n) satisfies a central 
and local limit theorem with mean 0.39089 n and variance 0.041565 n. 



Our particular strategy is rooted in our recent work on asymptotic enumeration of fc-noncrossing 
RNA structures [TT] and a paper of Bender [5] who showed how such central limit theorems arise 
in case of singularities that arc poles. In order to put our results into context let us provide some 
background on central and local limit theorems. Suppose we are given a set A n (of size a n ). For 
instance let A n be the set of subsets of {1, . . . , n}. Suppose further we are given A n ^ (of size a n ^. 
k G N representing a disjoint set partition of A n . For instance let A n ^k be the number of subsets 
with exactly k elements. Consider the random variable having the probability distribution 
P(£n = k) = a n> k/a n , then the corresponding probability generating function is given by 



P(£ n - k)w h = 2^ — w K = 

fc>n fc>n a ™ 2^,k>o a n,ki- 
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Let ip n (w) = X)/c>o a n,kW k , then ^" fej is the probability generating function of and 



f(z, w) = fn{w)z n = X! E a 



n,kW k Z T 



n>0 n>0k>0 



is called the bivariate generating function. For instance, in our example we have P(£ n ) = ipr an d 
the resulting bivariate generating function is 

x x ' ^w"z n — 



M EE 

n>0 fc<7i 



kj l-z(l + u>) 



The key idea consists in considering f(z,w) as being parameterized by w and to study the change 
of its singularity in an e-disc centered at w = 1. Indeed the moment generating function is given 

by 



fc>0 



¥>„(!) [*"]/(*, 1) 



and ^^Mw^jv ~ E{e lt ^ n ) is the characteristic function of £ n . This shows that the coefficients 
of f(z,w) control the distribution, which can, for large n, be obtained via singularity analysis. 
The resulting analysis can be amazingly simple. Let us showcase this in the case of the binomial 
distribution. Here we have the bivariate generating function X)«>o S 



k „n 



k<n \k) W Z — 1-z(1+wj) ' 

eq. CEO). The simple pole r(s) of f{z,e s ) is j^. Observe that ^-£2 - (^7§j) W holds for s 
uniformly in a neighborhood of 0, and Taylor expansion shows 

uniformly for t for any arbitrary finite interval. It remains to apply the Levy-Cramer theorem 
(Theorem 0]) to the normalized characteristic function of the random variable -^7=-, which yields 

the asymptotic normality of rj n . Thus (^) is asymptotically normal distributed with mean ^ and 
variance j. As it turns out we will have to work a bit harder to prove our main result. The 
complication is due to the fact that the generating function for 3-noncrossing RNA structures 
is much more complex (and fascinating) than the bivariate function of eq. {TTT]) which has a 
simple pole as dominant singularity. For instance, the singularity of the generating function for 
3-noncrossing RNA structures is not a pole but of algebraic-logarithmic type [16] . 



Our two main results, Theorem [5] in Section [4] and Theorem [6] in Section [5] shed light on the 
distribution of 3-noncrossing RNA structures from a global and local perspective. A central limit 
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theorem represents the global perspective on the limiting distribution of some random variable X n : 



lim 



X n — fj, n 



< x 



1 



<2n 



~dt 



Bender observed in [2] that a central limit theorem combined with certain smoothness conditions 
on the coefficients a n % implies a local limit theorem which considers the difference between P(x < 
Xn ~^ n < x + 1) and f* +1 e~^dt as n tends to infinity. To be precise, X n satisfies a local 
limit theorem on some set S C K if and only if 



lim sup 



A" 



1 



2tt 



= 



holds and we say X n satisfies a local limit theorem for some S = {x £ M | x = o(yfn)}. Why 
is the smoothness of the a n ^ so important? Suppose a n ,k = (^) + (— l) fe 2(^), then it follows in 
analogy to our above argument that a central limit theorem with mean jn and variance \n holds. 
However, r\ n does not satisfy a local limit theorem, since 



Vn 



1 



'In 



x)- 



1 



and for S = {^\n =1,2...}, we have ±V^P(?7 n = \ri) - -j= 



2tt 



flips between -(") and 3(") 



0, the key point being that 



All results of this paper hold for 2-noncrossing RNA structures, i.e. RNA secondary structures. 
This is a consequence of an analogous analysis of their respective bivariate generating function. In 
this case, however, no singular expansion is necessary as the generating function itself can be used. 
They also give rise to put the asymptotic results on RNA secondary structures of [5] on a new 
level. We can pass from computing exponential growth rates to computing distributions for RNA 
secondary structures with specific properties. To be precise we have for RNA secondary structures 



Theorem. Let S^n, h) denote the number of RNA secondary structures with exactly h arcs. Then 
the random variable Y n having distribution ¥(Y n = h) = S- 2 (n, h)/S2(n) satisfies a central and local 
limit theorem with mean 0.27639 n and variance 0.04472 n. 



In particular the theorem predicts a sharp concentration of the number of RNA secondary struc- 
tures with 55.278% unpaired bases which agrees with the statistics of RNA secondary structures 
obtained by folding algorithms [24l [HI [22l [TH [T3 [15]. Let us finally remark that much more 
holds: due to the determinant formula for fc-noncrossing matchings and the functional identity 
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FIGURE 4. Central limit theorem of 2-noncrossing and 3-noncrossing RNA structures: 
both random variables are normalized to S^n, /i)/S2(w)and S^n, h)/Sa(n), respectively. 
In case of n = 100, for 2-noncrossing RNA structures we have a mean of 0.276393 n = 
27.6393 and variance 0.044721 n = 4.4721 (left curve), while for 3-noncrossing RNA 
structures mean 0.39089 n = 39.089 and variance 0.041565 n = 4.1565 (right curve). The 
red dots and magenta dots represent the values S' 2 (?i, h)/S2(n) and S' 3 (n,h)/S3(n), re- 
spectively. 

of Lemma [TJ Section [3] our results can be generalized to fc-noncrossing RNA structures, where fc 
is arbitrary. Why this is of interest can be seen in Fig. [4] For higher fc the mean of the central 
limit theorems for fc-noncrossing RNA structures will shift towards the maximum combinatorially 
possible number of arcs. We speculate that each increase in fc will basically cut the distance to the 
maximum arc number in half. This is work in progress. 

The paper is structured as follows: In Section [2] we provide some background on fc-noncrossing 
RNA structures and all generating functions involved. In Section [3] we give a functional equation 
for the bivariate generating function of S' 3 (n, h) via 3-noncrossing matchings proved in [llj . We 
have included its proof in the appendix in order to keep the paper self-contained. This functional 
identity plays a key role in proving the central limit and local limit theorem in Section U and 
Section [5l respectively. The central limit theorem is proved by analyzing the singular expansion of 
analytic function of power series X)n>o Eft<a ^(n, h)w z n and using transfer theorems [5l l6l [16] 
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and to prove the local limit theorem, we use a theorem of Hwang [5] and build on our proof of the 
central limit theorem. 



2. RNA STRUCTURES 

Let us begin by illustrating the concept of RNA structures. Suppose we are given the primary 
sequence 

AACCAUGUGGUACUUGAUGGCGAC . 

Structures are combinatorial graphs over the labels of the nucleotides of the primary sequence. 
These graphs can be represented in several ways. In Figure [5] we represent a 3-noncrossing RNA 
structure with loop-loop interactions in two ways: first we display the structure as a planar graph 
and secondly as a diagram, where the bonds are drawn as arcs in the positive half-plane. 

In the following we will consider structures as diagram representations of digraphs. A digraph D n 
is a pair of sets VD n ,En n , where Vd„ = {1, .. and En n C | 1 < i < j < n}. Vd„ and 

Eu n are called vertex and arc set, respectively. A A:-noncrossing digraph is a digraph in which all 
vertices have degree < 1 and which does not contain a fc-set of arcs that are mutually intersecting, 
i.e. 

(2.1) (^ri)Jri)) (^2? ^2)1 • ' ' 3 (^rfcjjVfc)) Vi ^ < ' ' ' < Vfc ^ jri ^ 3ti ^ ' ' ' ^ J^k ' 

We will represent digraphs as a diagrams (Figure [5]) by representing the vertices as integers on a 
line and connecting any two adjacent vertices by an arc in the upper-half plane. The direction of 
the arcs is implicit in the linear ordering of the vertices and accordingly omitted. 

Definition 1. An RNA structure (of pseudo-knot type k — 2), is a digraph in which all 

vertices have degree < 1, that does not contain a fc-set of mutually intersecting arcs and Tares, 
i.e. arcs of the form + 1), respectively. We denote the number of RNA structures by Sfe(n) 
and the number of RNA structures with exactly £ isolated vertices and with h arcs by Sfe(n, £) and 
S' k (n, h), respectively. Note that S' fe (n, h) = Sk(n, n — 2h). 
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FIGURE 5. A 3-noncrossing RNA structure, as a planar graph (top) and as a diagram (bottom) 



Let fk(n,£) denote the number of fc-noncrossing digraphs with I isolated points. We have shown 
in [TU] that 



(2.2) 



fk(n,£) = ( - e )f k (n-l,0) 



T 

(2.3) det[J i _ J (2aj) - 7 i+J (2x)]|^ 1 = ]T f k (n, 0) • — 



(2.4) e -dct[/ l _,(2x)-/ l+J (2 a; )]|^ 1 = (^^)(^/ fc (n,0): 

l>0 ' n>l 

In particular we obtain for k = 2 and k = 3 



77 \ / 77 

(2.5) f 2 (n,£)= ( JC M)/2 and f 3 (n,£) 



EEam 



n>l kf=0 



~ — 2~~ - o rl 



where C m denotes the m-th Catalan number. The derivation of the generating function of k- 
noncrossing RNA structures, given in Theorem [T] below uses advanced methods and novel con- 
structions of enumerative combinatorics due to Chen et. al. [H [7] and Stanley's mapping between 
matchings and oscillating tableaux i.e. families of Young diagrams in which any two consecutive 
shapes differ by exactly one square. The enumeration is obtained using the reflection principle due 
to Gessel and Zeilberger [7] and Lindstrom [13] combined with an inclusion-exclusion argument 
in order to eliminate the arcs of length 1. In [lOj generalizations to restricted (i.e. where arcs of 
the form (i,i + 2) are excluded) and circular RNA structures are given. The following theorem 
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provides all data on numbers of fc-noncrossing RNA structures with h arcs and the numbers of all 
fc-noncrossing RNA structures. 

Theorem 1. [10] Let k G N, k > 2, let C m denote the m-th Catalan number and fk(n,£) be the 
number of k-noncrossing digraphs over n vertices with exactly £ isolated vertices. Then the number 
of RNA structures with £ isolated vertices, Sk(n,£), is given by 



b=o v 7 

where (n — 2b, £) is given by the generating function in eq. (|2.3[) . Furthermore the number of 
k-noncrossing RNA structures, Sfc(ri) is 



where {X)"=o fk(n — 2b,£)} is given by the generating function in eq. (|2.4p . 

In principle, Theorem [T] contains all information about the numbers of /c-noncrossing RNA struc- 
tures. However, due to the inclusion-exclusion structure of its coefficients it is however difficult to 
interpret and to express their behavior for large n. Subsequent asymptotic analysis |11| produced 
the following simple formula 

Theorem 2. [11] The number of 3-noncrossing RNA structures is asymptotically given by 



We have shown in the introduction that the bivariate generating function is the key to prove the 
central and local limit theorems. The following lemma, whose proof is given in the appendix, 
rewrites this bivariate generating function as a composition of two "simple" functions. This is 
crucial for the singularity analysis insofar as we can use a phenomenon known as persistence of the 
singularity of the "outer" function (the supercritical case) [5] . It basically means that the type of 
the singularity is determined by the generating function of /c-noncrossing matchings. 



(2.6) 




(2.7) 





3. A FUNCTIONAL EQUATION 
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Lemma 1. Let x be an indeterminant over R and w G K a parameter. Let pk{w) denote the 
radius of convergence of the power series Eh>dEkh/2 S' fc (n, h)w 2h ]x n . Then for \x\ < pk{w) 

(3.1) v v sk„,/.XV = ^_^ EA(2 „,o) (j^^) 

n>0h<ri/2 n>0 x ' 

holds. In particular we have for w = 1 

1 / \ 2ri 

(3.2) E = E 0) (^n) 
/or zgC with \z\ < /%(!). 



To keep the paper selfcontained we give the proof of Lemma Q] in the Appendix. While <|3 . 1 1) 
can only be proved on the level of formal power-series for real variables, complex analysis i.e. the 
interpretation of these generating functions as analytic functions allows to extend the equality to 
arbitrary complex variables. 

Lemma 2. Suppose e > 0, k £ N, k > 2 andw — ei , where \s\ < e and (p n ,k(s) — Tl,h<n/2 ^'k( n ' h)e hs . 
Let Pk(s) G R + denote the radius of convergence of^2 n>0 (p nt k(s)z n parameterized by s. Then we 
have 
(3.3) 

VveC; |s| < e, \z\ < Pk {s); ^ <fin,k(s)z n = 2 j ]T / fc (2n, 0) 

n>0 7i>0 

Furthermore X) n >o L Pn,3{s)z n has an analytic continuation, ^{z,s). For e sufficiently small and 
\s\ < e, 23(2:, s) /ias exactly 6 singularities, 4 of which have distinct moduli. 



( e2 z \ 
\e s z 2 — z + 1 / 



Proof. We first prove eq. (|3.3p . For this purpose we observe that 
(3.4) V \s\ < e, \z\ < Pk (s) G(z, s) = 1 E /*( 2 »- °) » 

n>0 



considered as a power series in e^ s is analytic in a neighborhood of s = 0, since G(z, 0) is analytic 
for |z| < pfc(0). In addition, we can interpret J2 n >o L Pn,k(s)z n as a power scries in e^ 8 : 



(3-5) E E 



n>0 \/i< 



r s; (n, ^) e ' is ] z n = J2 ( E s *(n, h >" 1 ( eS )' 1 = E M*) ( ei 

n/2 / h>0 \n>2h J h>0 



21 1 
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Therefore G(z 1 s) and the power series X) n >o L Pn,k{s)z n are analytic in the indeterminant e^ s in 
an e-disc centered at 0. Lemma [T] implies that for s g] — e, e[ the analytic functions G(z,s) and 
Xm>o L Pn.k{s)z n are equal. Since any two functions that are analytic at and that coincide on the 
interval ] — e, e[ are identical, we obtain 

(3.6) V \s\ < e, \z\ < Pk (s) G(z, s) = £ cp n ,k( s >" ■ 

Claim 1. Suppose |s| < e. Then ^ re >o W3(s)z" has an analytic continuation, 23(2, s), which has 
exactly 6 singularities 4 of which have distinct moduli. 

In order to prove Claim 1 we observe that the power series J2 n >o f^i^ 71 ' 0)2/™ has the analytic 
continuation vE'(y) (obtained by MAPLE sumtools) given by 

-(i-ie^ip^C-^) 

(3-7) = — -4 ^- , 

16 J/2 

where P™{x) denotes the Legendre Polynomial of the first kind with the parameters v = | and 
m = — 1. According to cq. (|3.6[) we have 

, » v 2n 

(3-8) E = ^TTT E ^ 2 "' °) (^Stt) 

n>0 n>0 V 7 

which implies that X) n >o l Pn,z{s)z n has the analytic continuation 
(3.9) V| S |<6, E 3 (z,s)= - 5 1 — 



e s z 2 — z + 1 I e s z 2 — z + 1 J 

In particular for s = 0, 23(2, 0) is the analytic continuation of the power series X)n>o 53(71)2". We 
proceed by showing that 23(2, s) has exactly 6 singularities and 4 of them have different moduli in 
C parameterized by s. Two singularities are given by the roots of e s z 2 — z + 1 are Ci( s ) = 1 ~X^7 4e - 
and ( 2 (s) = 1+ ^I 4eS ■ Observe that |Ci(0)| = | <" 2 (0) | = 1 and polynomial e s z 2 -z + 1 depends 
continuously on ei , therefore Ci( s ) an d £2(3) could potentially have equal modulus for |s| < e. The 
remaining 4 singularities are induced by the the unique dominant singularity a\ = of analytic 
function ^>(y). The function ^(y) has three singularities, two of them ol\ = jg and a.2 — +00 are 

branch points and the other 0:3 = is a removable singularity. The function g(z) = ( e e JiL^+i j 
with g(0) = has a radius of convergence of 1 as s tends to 0. Therefore the singularity type only 
depends on ^>(y) (this is the supercritical case in [5]). The singularity a.\ = gives rise to the 
equations 

= e/z 2 - (l + 4e^ s )z + l and = e s z 2 + (4,ei s - l)z + 1 



2 
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and setting p+(s) = 1 + 4e? s , fi-(s) = 1 — Ae^ s and 9{s) = y/ I2e s + 8ei s + 1 its roots arc given 

by 

respectively. Observe that for \s\ < e, e 5 is in a neighborhood of 1 over C, hence 9(s) ^ 0. That 
leads to 4 distinct roots ^3(3), £4(5), ^5 (s), ^(s) over |s| < e, all of them have distinct moduli for s 
being a sufficiently small neighborhood of 0. Indeed, for s = we have 4 distinct real valued roots 

Ca(0) = g , Ct(0) = ^ ' &( ) = 2 ' ^ ^ = 2 

and the polynomials e s z 2 — (1 + 4e3 s )z + 1, e s z 2 + (4e3 s — l)z + 1 and e s z 2 — z + 1 depend 
continuously on the parameter ez 8 , whence Claim 1 and the lemma follows. □ 



4. The central limit theorem 



In this section we prove a central limit theorem for the numbers of 3-noncrossing RNA structures 
with h arcs. We will analyze for fixed but arbitrary n the distribution of S' 3 (n,h). Let us first 
prepare some methods and results used in the proof of Theorem[5] [z n ] f(z) denotes the coefficient 
of z n in the power series expansion of f(z) around 0. The scaling property of Taylor coefficients 

(4.1) V 7 £<C\0; [z n ]f(z)= 7 n [z n ]f(~) , 

1 

shows that w.l.o.g. any singularity analysis can be reduced to the case where 1 is the dominant 
singularity. We will be interested in the behavior of an analytic function "locally" , i.e. around a 
certain singularity p. For this purpose we use the notation 

(4.2) f(z) = O (g(z)) as z — > p -t=^ f(z)/g(z) is bounded as z — > p 

and if we write f(z) = 0(g(z)) it is implicitly assumed that z tends to a (unique) singularity. 
Given two numbers <f>, R, where R > \p\ > and < cj) < ^ and p £ C the open domain A p (0, R) 
is defined as 

(4.3) A p (cf>, R) = {z\ \z\ <R,z^p, |Arg(z - p)\ > 0} 

A domain is a A p -domain if it is of the form A p (c/>, R) for some R and <fi. A function is A p -analytic 
if it is analytic in some A p -domain. We use U(a,r) = {z £ C | \z — a\ < r} to denote the open 
neighborhood of a in C. Via the following theorem we can extract the coefficients of analytic 
functions provided these functions satisfy certain "local" properties. 
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Theorem 3. [5] Let r £ Z>o and f(z,e s ) be a A p ( s y analytic function parameterized by s, which 
satisfies in the intersection of a neighborhood of p(s) with its A p r s \- domain 

(4.4) f(z,e s ) = b (s)+b 1 ( s )(z-p( s )) + A( s ) (p(s) - z) r ln ( ~rj^ ) + R{z,a) 

where A(s),ba{s),b\{s) are analytic in \s\ < e and \R(z,s)\ < c\p(s) — z\ for some absolute constant 
ceC. That is we have f(z,e s ) = 0((p(s) — z) r ln( p ^_ z )) with uniform error bound as s in a 
neighborhood o/O. Then we have 

(4.5) [z n ]f{z, e s ) = A(s) (-l) r - ^— ( 1 - 0(-)) for some A(s) £ C , 

n(n — 1) . . . (n — r) \ n ) 

where the error term is again uniform for s from a neighborhood of origin, i.e. R(s) < c \ s\, where 
c> 0. 



Remark. The equivalence between eq. (|4.4p and f(z,e s ) = 0((p(s) — z) r m ( p ^*_ z )) for r E 1>q 
can be seen as follows: by definition of f(z, e s ) = 0((p(s) — z) r ln( p ^_ z )) there exist A(z, s) and 
B(z, s), such that f(z, e s ) — B(z, s) + A(z, s)(p(s) — z) r hi( p ( g )_ z ); where A(z, s) and B(z, s) are 
analytic in a neighborhood of p(s). Taylor expansion of A(z, s) and B(z, s) at z = p(s) produces 

f(z, s) = B(z, s) + A(z, s)(p(s) - z) r In 

= b Q (s) + h(s)(z - p(s)) + ■■■ + (a (s) +ai(z- p(s)) + ■■■) (p(s) - z) r \n 
= b (s) + bi{s)(z - p(s)) + a {s)(p{s) - z) r In ^ 1 - z ) + R ^ ^ 

where R(z,s) = 0((p(s) - z)^ In (^) ). For r e Z> , gg^- = 0(\p(s) - ,|nn( R ^)) 
is bounded by an absolute constant as z tends to p(s). That implies the error bound is uniform. 

The next Theorem is a classic result on limit distributions which allows us to prove our main result 
via characteristic functions i.e. explicitly by showing limn—^ ip n (t) = ip(t) for any t G (—00,00). 

Theorem 4. (Levy-Cramer) Let {£ n } be a sequence of random variables and let {(p n (x)} and 
{F n (x)} be the corresponding sequences of characteristic and distribution functions. If there exists 
a function (p(t), such thatlmin^oc ip n (t) = <p(t) uniformly over an arbitrary finite interval enclosing 
the origin, then there exists a random variable £ with distribution function F(x) such that 

F„ (t) =^ F{x) 

uniformly over any finite or infinite interval of continuity of F{x). 
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We now consider the random variable X n having the distribution F(X n = h) = S' 3 (n, h)/S3(n), 
where h = 0, 1, . . . [_§ J ■ The key point in the proof of Theorem [5] is to compute the coefficients of 
the bivariate generating function whose variable, s is considered as a parameter. Intuitively the 
particular distribution is a result of how the singularity shifts as a function of this parameter. As 
a result the proof is somewhat "non-probabilistic" and has two distinct parts: (a) the analytic 
combinatorics of the bivariate generating function and (b) the computation of the characteristic 
function with subsequent application of the Levy-Cramer Theorem. 

Theorem 5. The random variable ~ 2 ^ n has asymptotically normal distribution with parameter 
(0,1), i.e. 

(4.6) lim P ( < x) = -L f e-^ d t 

n->-oo y VCT 2 ?i / V27T J-oo 

and /x, a 2 are given by 

4.7) fj, = | 42 ;_ = 0.39089 and a 1 = u 2 ^= — = 0.041565 . 

2 2 V ix 2 



Proof. We set w = e2 s and (p n ,3(s) = J2h<n/2 %( n ' h)e . Since 

(4.8) J2 <p«*(*)* n = E ( E s 3(^' ^ I *" . 

n>0 n>0 \/t<n/2 / 

we can consider the double generating function En>o EKn/2^(' ! '' i )''' 2 '' z ° as a P owcr scries in 
the complex indetcrminant z, parameterized by s. 
Claim 1. 



(4.9) = O [{l - 16z) 4 In ( jTTiei)) 



holds uniformly for Vz G A^^, JJ) n ^(^; e ); 



^(z) is A j_ (0, i?)-analytic and has the singular expansion (1 — 16z) 4 In I i_\q z ) m the intersection 
of U(jjr, e) with the A_^— domain, where A r ((f>, R) = {z \\z\ < R, z r, |Arg(z — r)| > </>} for some 
i? > r. First A j_ (</>, i?)-analyticity of the function (1 — 16z) 4 ln (^ j^^ ) is obvious. Wc proceed 

by proving that (1 — 16z) 4 ln f TZig; ) is the singular expansion of ^(z). The above mentioned 
scaling property of Taylor coefficients allows us to consider the power series J2 n >o f^i^ 71 ^ ®)(fe) n 
over the A-domain Ai(</>, R) for some R > 1. Using the notation of falling factorials (?i — 1)4 = 



16 
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(n — l)(n — 2)(n — 3)(n — 4) we observe 



f 3 (2n,0) = C n+2 C n -C 2 n+1 = 



1 



12(n- l) 4 (2n+ 1) (2n 



(n - 1) 4 (n + 3)(n + l) 2 {n + 2) 2 V n 
With this expression for /3(2n,0) we arrive at the formal identity 



^16-"/ 3 (2n,0)z" = 0(J2 



n>5 



n>5 



16 _n 



12(n- l) 4 (2n + 1) (2n 
n 



(n-l) 4 (n + 3)(n + l) 2 (n + 2) 5 



4! 1 1 



(n — 1)4 7T 71 



(n — 1 4 7r n 



n>5 



where f(z) = 0(g(z)) denotes that the limit f(z)/g(z) is bounded for z — * 1, eq. (|4.2p . It is clear 
that the error bound below 



E 

n>5 

E 

n>5 L 



16- 



16" 



1 



12(n- l) 4 (2n + l) (2n 



(n - 1)4 (n + 3)(n + l) 2 (?i + 2) 2 V, n / (n - 1) 4 tt n 



4! 1 1 



12(n- l) 4 (2ra+ 1) (2n 



(n - 1)4 (n + 3)(n + l) 2 (n + 2) 2 \ n / (n - 1) 4 7r n 



4! 1 1 



< K 



holds uniformly for z in Ai(<f),R) n J7(l,e) and some absolute k < 0.0784. Therefore we can 
conclude 



(4.10) 



£ 16-'7 3 (2n,0)z" = G>(£ - A-il^) . 
- - (n — 1)4 tt n 



n>5 n>5 

We proceed by interpreting the power scries on the rhs, observing 



(4.11) 



Vn>5 



z n ] (1-z) 4 In. 



1 



4! 



1 



1 — z J (n — 1) ... (n — 4) n 



whence ^(1 — z) 4 In y^-j^ is the unique analytic continuation of J2n>5 (n-i) 4 h n 2 ™' Using the 
scaling property of Taylor coefficients [z n ]f(z) = J n [z n ]f{—) 



(4.12) V(z) = (l-16z) 4 ln 



1 



1 - 16z 



holds uniformly for V z € A_i_ (0, i?) n ^(jg > e ) 



16 ' 



Therefore we have proved that (1 — 16z) 4 ln( 1 _ L 16z ) is the singular expansion of ^f(z) at 
whence Claim 1. Our next step consists in verifying that when passing from *&(z) to the bivariate 
generating function &(z, s) = ^{{ w 2 z fl z+l ) 2 ), then there exists a singular expansion of the form 
O ^(1 — ^j^j) 4 ln( 1 1 5 - )j , parameterized in s. 

Claim 2. Let < e < 1, then for any |s| < e and 2 g A p3 / s -\((f),R), we have ^(z, s) = 
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^(1 — pa * s j ) 4 ln( 1 * an d the error bound is uniform for s in a neighborhood of 0. 

To prove the claim we first observe that Claim 1 implies 



(4.13) 



#0) = k(1 - 16z) 4 ln 



1 



1 - 16z 



i?(z) 



for some absolute constant k and is the uniform error bound for z £ A^_(cf>,R) n [/(yg,e). 
I.e. For z £ A_i_ (</>, R)flU(jg, e), there exists some absolute constant c, such that |i?(z)| < c- 1 1 — 16z\ 
holds. According to Lemma [5] we have 



E 3 (z,s) 




+ R(z,e s ) 



We expand I 1 — 16( 



tion of 



+t) 



4 

2 ^ In 



1-16(- 



~n — 

i2 z 



around z = p 3 {s), where pz(s) is the Sein- 



ed 2 — z+1 



j of minimal modulus. Lemma [2] implies that 



Pais) = 



4ei s + 1 - Vtte 



2e s 



is the unique dominant singularity. As a function in s we have p' 3 (0) 



21 ^ 0. The term 



and to 



interval between 2 In i 

6 



1 in ps(s) produces two branching points parameterized by s. i.e. w = e^ s = — g 
— 5, or equivalently s = 2 In i + 27ri and s = 2 In g + 2ni, respectively. The 
2ni and 2 In i + 2-7TJ divides the complex plane of s into two analytic 
branches. For any < e < min{|2 In | + 2m\, \2 In ^ + 27ri|} = 6.4343, the region |s| < e is disjoint 

to the interval [(2 In |, 2tt), (2 In |, 27r)]. Therefore P3(s) is analytic for |s| < e. We next consider 

1 s 

q(z, s) = 1 — i6( e s z e 2_^ +1 ) 2 as a function of z and compute the Taylor expansion at P3(s). 



q(z, s) = a(p 3 {s) - z) + 0{z - p 3 (s)) 2 



18 
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and setting a = - 

5— v^l 

(a(p 3 (s) -z) + 0{z - p 3 (s)) 2 ) 4 In a{p3{s) _ z) ^ 0{z _ p3{s)y2 



e s z 2 — z + 1 



q(z,s) In 



q(z,s) 



e s (z - p 3 (.s)) 2 + (2p 3 (s)e s - l)(z - p 3 (s)) - 3p 3 (s)V + p 3 (s) + 1 

([a + Q(z - p 3 (s))](P3(s) - ^) 4 ln [^.p,; 1 ,,)],,,,,),,) ) 
O(z- A3 ( s ))-3 (03 ( S ) 2 +p 3 ( S ) + l 

0( P3 ( s )-z) 4 h/ 



p 3 (s) - z 

According to Theorem[3]for r = 4, we obtain the error term in the expansion of eSz <y_ z+l q(z, s) 4 In ^ * s ^ 
is uniform for s in a neighborhood of 0. We observe that the resulting error bound for S 3 (z, s) is 
the sum R(z, e s ) + R\{z, e s ), where 

\R(z, e*)| <c|l-16 ( s J^ Z z + 1 ] I = 0(p 3 (s) z) . 

Therefore the error bound for the expansion of bivariate S 3 (z, s) is uniform and Claim 2 is proved. 
We proceed by using the scaling property of Taylor coefficients \z n ]f(z) = 7"[ z ™]/(^) an d apply 
Thcorcm[3] Via Thcorcm|3]we obtain the key information about the coefficients of S 3 (z, s) which 
allows us to substitute <p nt3 (jf-) in eq. (|4.17[) below: 

(4.14) [z n ] E 3 (z, s) = K(s) ^(s)- 1 )" (l - 0(-)) for some K(s) G C , 

n(n — 1) . . . [n — 4) v ' \ n J 

where the error term is again uniform for s from a neighborhood of origin. 

Suppose we are given the random variable (r.v.) with mean [i n and variance a 2 . We consider 
the rescaled r.v. rj n = (£„ — ^, n )a~ l and the characteristic function of rj n : 

(4.15) /„„ (t) = E[e*"»] = E^^je"^* . 
In particular, for £„ = A"„ we obtain, substituting the term E[e z *' ,T '] 

(4-16) fxM - (± S ~^e^ . 

Since <p n , 3 ( s ) = E/ i < n /2 S 3( fi .' l ) e '' i . wc can interpret S 3 (n) = Z) A <„/ 2 S 3 (n, /i) as <£ n , 3 (0) and 

/2 C 



</'„ j3 (^-) = 53/i<n/2 S 3 (n, h)e h "" , respectively. Therefore we have 



(4-17) f Xn (t)=* iPn (!i)e- i % t 
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For \s\ < e, eq. yields <p(s) = [z n ]E 3 (z,s) ~ K(s) „ (n _ 1) 4 '. ( „_ 4) 

error term and we accordingly obtain 



with uniform 



(4.18) 



/*„(*) 



K(0) 



P3(0) 



where the error term is uniform for t from any bounded interval. Taking the logarithm we obtain 
(4.19) 



In fx n (t) ~ In ; : CT ", - — n In 



X(0) 



Ps(0) 



z — c 



Expanding g(s) = In in its Taylor series at s = 0, (note that g(0) = holds) yields 



(4.20) In 
and therefore 



K(^) 

(4.21) ln/„(t)~ln \ -n 



Psj^) = p 3 (0) ft 
p 3 (0) p 3 (0)a„ 

p 3 (0) it 



Ps'(0) /P' 3 (0) 



p 3 (o) Vp 3 (o) 



^ 2 

2(T 



S3 + 0((-J ) 



P 3 (0) er„ 



P 3 '(0) /V 3 (0) 



p 3 (o) Wo) 



Claim 2 implies 5 3 (z, s) = O KP3(s) — z) 4 In p3 (^_ z ^ is analytic in s where s is contained in a 
disc of radius e around 0. Hence S 3 (z,s) is in particular continuous in s for |s| < e and we can 
conclude from eq. (|4.14[) for fixed t &] — oo, oof 



(4.22) 



lim ( InK(-) -ln/i(0) ] = 



In view of eq. (|4.2ip we introduce 



and eq. (|4.2ip becomes 
(4.23) 



p 3 (0)' 



ln/x„(<) 



p 3 (o) 



0( 



p 3 (o) 



with uniform error term for t from any bounded interval. This is equivalent to linin^oo fx n (t) 
exp(— y) with uniform error term. The Levy-Cramer Theorem (Theorem U) implies now eq. 



and it remains to compute the values for /x and a which are given by 

p'M 



(4.24) 
(4.25) 



P 3 (0) 



2 42 - 0.39 89 



I-IV21 



p 3 (o) 



441 ■= 0.041565 



5-V21 
2 



2(1 
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whence eq. (|4.7[) and the proof of Theorem [5] is complete. □ 



5. The local limit theorem 



In this section we complement the central limit theorem presented in the previous section 

lim P I Xn ~ g» < x] = __L 



n->oo V o n I J2TT J_ 



oo 



by considering a "local" perspective on the limiting distribution of X n , For the local limit theorem 
we analyze the difference between ¥(x < jY " g ~ M ' 1 < x + 1) and —1= e~^dt as n tends to 
infinity. X n satisfies a local limit theorem on some set S C K if and only if 



(5.1) lim sup 



X n ~ (J"n \ 







One key condition formulated in cq. (|5.2[) of Theorem [6] below for proving a local limit theorem is 
given by 

^±~ eX v(M{s)P n +N(s)) , 

<£n(0) 

where M(s) is differentiable and N(s) is continuous in some e-disc centered at 0. In view of 
eq. (|4.17p and eq. (|4. 18|) in the proof of the central limit theorem, this condition alone implies the 
central limit theorem. In other words, the local limit theorem implies the central limit theorem. We 
have shown in the introduction that a central limit theorem does not imply a local limit theorem. 
Bender observed in [2] that the central limit theorem combined with certain smoothness conditions 
docs imply the local limit theorem. Accordingly, in order to prove the local limit theorem for 3- 
noncrossing RNA structures with h arcs our strategy will consist in verifying such smoothness 
conditions [5]. 

Theorem 6. Let <p n {s) = J^ fc a n: kW k and w = e s . Suppose 
(5.2) ^ ~ exp(M( S )/3„ + N(s)) 

holds uniformly for \s\ < r, s G C and t > 0, where the following conditions are satisfied 

(i) M(s) is differentiable and N(s) is continuous in \s\ < e and furthermore M{s) and N{s) are 
independent of n. 

(ii) f3 n is independent oft, f3„ — ► oo and M"(0) > 0; 
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(iii) there exist constant 5 and c = c(<5, r) > ; where < S < r such that 

<Pn{r + it) 



(5.3) 



0(exp(-c/J„)) 



<Pn(r) 

holds uniformly for — r < r < r and <5 < |t| < 7r as n tends to infinity. 

Then random variable X n having distribution ¥(X n = k) = a n ^ja n with mean M'(0)/3 n and 
variance M"(0)j3„ satisfies a local limit theorem on the real set S = {x \ x = o(y / /5^)} i.e. 



(5.4) 



lim sup 



X n — fj, r , 



1 



= 



With the help of Theorem [6] we can now prove the local limit theorem for 3-noncrossing RNA 
structures with h arcs. 

Theorem 7. Let S' 3 (n,h) be the number of 3-noncrossing RNA structures with exactly h arcs. Let 
X n be the r.v. having the distribution 

(5.5) Vh = 0,l,...£\, F(X n = h) = ^^ 

2 S 3 (n) 

Then we have for set S = {x \ x = o(^/n)} 

(5.6) lim sup 
where pi = 0.39089 and a 2 = 0.041565. 




Proof. We will show that { ^f^F } satisfies the conditions for Theorem [5] For \s\ < e, where e is 
sufficiently small but fixed. The crucial equation implying the conditions of Theorem [5] is eq. (|4.14|) 
of the proof of Theorem [5] 

<^, 3 (s) = K(s)- -Ms)- 1 )* 1 (l O(-)) K(s) e C , 

n(n — 1) . . . (n — 4) \ n J 

holds uniformly for |s| < e. Therefore we have 

Pn,3(0) K(0)\ P3 (s)J \ y n'J P V W*)/ \K(0))J 
uniformly for |s| < e. We set 

(5.8) (3 n = n, M(s) =ln(^) and N(s) = In . 
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By construction t is independent of n and clearly n — > oo and M(s) is differentiable and N(s) 
is continuous for all s such that \s\ < e. In addition M"(0) is analytic for \s\ < e and we have 

1 _ 94 /qT 

M"(0) = /J, 2 - 5 4 _% = 0.041565 > 0. Let S = e and -e < r < e, we obtain 



^n,3(^) 



cxp n 



In 



P3(r) 
Pz{r + it) 



In 



K(r + it) 
W) 



uniformly for — e < r < e and e < |t| < n. Since yields a constant factor and K(s) is continuous 



for \s\ < e, it suffices to analyze In 



Pa(r) 
P3(r+it) 



K(0) 

We observe p 3 (s) 

P 3 (0)' 



l+4e2 



V 



12e s 



complex s where \s\ < e. The singularities of ln( ^| s j ) correspond to the zeros of 12e s 



^ for any 
- 8e* + 1 = 



(2es + l)(6e2 + 1), that is = — i or — g. Observe that for |s| < e, | e S | is close to 1. Therefore 



In 



is analytic for any e < |t| < n and r g] — e, e[ and we can conclude 



O 



exp(?i • In 



P3(it) 
Ps(P) 



= O I exp Re u ■ In 



p 3 (r + it) 
Ps(f) 



uniformly for — e < r < e and e < \t\ < 7r. Taylor expansion of m ( P3 pJ(^t ) a ^ shows (see 
cq. (|4.20| ). that the dominant real part of ln( p^iy) ) ^ s gi ven by 



Ps(r) 



2! 



< for r e] - e, e[ 



Setting ci = 



P3 ( r ) 




2 

7T 


P3<>) 




2! 



Ps(r) 


/PaM^ 2 " 


5 2 


P3(r) 


WO J 


2! 



</V3(V + it) 



(9(exp(— c ■ n)) 



<^n,3(r) 

for some < C2 < c < ci, uniformly for — e < r < e and e < |t| < it and Theorem [6] applies, whence 
Theorem [3 □ 



6. Appendix 



Proof of Lemma [TJ First we observe that for x,w E [ — 1, 1] the term w 2 x 2 — x + 1 is strictly 
positive. We set 

(6.1) F k (x,w) = J2 E %{n,h)w 2h x n 

n>0 h<n/2 
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and compute 

j I \ 2(h - j) 

n>0h<n/2j=0 v J 7 v v J ' 

n / 2 / \ / o- 



f»(x,») = EED- 1 )'(";1(,I ! i)''»-'i»i'"'*'' 



e e E(-')i"; J j( 2 ;"4j«^-A»)^ 



n>0 j<n/2 h=j 

n/2 , -\ / r» ■ 

n-j\ n-2j 



j>0 n>2j h=j 



e e E(-D j ( ; y j ; } ) M2(h - j), o )w ^ 



= »-i^E(»-i)'E^jA(2(*-i),o) (K 

j>0 J n>2j fi=j V V JJ/ V J; 

Wc shift summation indices n! = n — 2j and h' = h — j and derive for the rhs the following 
expression 



B-D'^ £ («' + D' £ ( 2( ;: ,) ap(» - m w ^'-" 



n/2 

E 

j>0 n'>0 ft=j 

E( - 1)) (^ E( „. +J)! ('' ,2 £'' , V" v ) A(2y ,„ ) ^ i *• 

j>0 J ' n'>0 h'=0 ^ 7 



The idea is now to interpret the term J2h'=o (2k') fk(2h' \ 0)w 2h ^ as a product of the two power 
«>oA(2n,0)^" 

\2n 



series and E„>o/*( 2jl ' ) — 



£|£A(2»,o><f^ = E £ {^A^o^k 



£>0 n>0 n'>02n+i=n 

'n'/2 



e e ; A( 2 »,o)»- . 

n'>0 I n=0 v 7 



We set ?7„' = {XX=0 (w) f k fth' , 0)w 2 ' l '|. By assumption we have |x| < pk{w) and we next 
derive, using the Laplace transformation and interchanging integration and summation 

(6.2) (»' +i)'^^r = /°° E Vw^Ve-tdt . 



ra'>0 u n'>0 
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Since |x| < pk(w) the above transformation is valid and using 

(-) E {e(^.^)^-E^E/.(^»)^ 

n'>0 {n=0 v 7 J £>0 n>0 v ; 

we accordingly obtain 

(6.4) £^£j£W«dt = /°°e-5:/ fc (2n I 0)^p^e-^. 

n'>0 n>0 ^ 

The next step is to substitute the term J2n'>o( n ' + J ) ' ^n. 7PT ' n e< l" ^.2p . whence consequently 



j>0 ''' ^° n>0 



j>0 J ra>0 

The summation over the index j is just an exponential function and we derive 

\2n 



f i 



n>0 

1 / \ 2 ™ 

' " 2^,2 _ , iu\2n, 



w 2 x 2 — X + 1 



((w 2 x 2 - x + l)ty n dt 



n>0 



We proceed by transforming the integral introducing u = (w x — x+l)t, i.e. dt = (w x — x+1) du 
and accordingly arrive at 

1 / \ 2n POO 

n>0 

= EA-( 2 ^°)7^if 2 2^ ^t ) 2 " 2 2 1 ^ ( 2 ") ! 

•'-^ (2n)! V urar — x + 1 J wrar — x + 1 

n>0 V / \ / 

27? 



MJ 2 X 2 — .T + 1 



n>0 

In particular for it) = 1 

(6.5) ]TS fc (n)x" = ^-L— ^)/ fc (2n,0) 



x — x + 1 ^-^ \ i - 1 + 1 

n>0 ri>0 



2rx 

X 



holds for any x £ M, satisfying |x| < and where pfe(l) is the radius of convergence of the 

power series J2n>o->k(n)z n over C, that is cq. (|6.5[) holds for x £] — Pfe(l), Pfc(l)[ ■ From complex 
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analysis we know that any two functions that are analytic at and coincide on an open interval 
which includes are identical. Therefore eq. (|6.5[) holds for z € C, \z\ < Pfc(l), and the proof of 
the lemma is complete. □ 
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